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FOREWORD 


This  technical  summary  report  was  prepared  by  the  Research  Department  of  the  Allison  Divi¬ 
sion  of  General  Motors  Corporation.  The  work  reported  was  accomplished  under  Contract 
Nonr-4  104  (00). 

The  program  was  sponsored  by  the  Advanced  Research  and  Project  Agency  through  the  Power 
Branch  of  the  Office  of  Naval  Research  under  the  direction  of  Dr.  J.  Huth  of  ARPA  and  Mr. 

J.  A.  Satkowski  of  ONR. 
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I.  INTRODUCTION 


This  third  Quarterly  Technical  Summary  Report  describes  the  progress  made  in  the  period 
1  November  through  3!  January.  During  this  period  the  MPD  section  was  incorporated  into 
the  closed  loop  device.  A  large  number  of  runs  with  pure  helium  were  made.  Also,  four  runs 
were  made  with  helium  seeded  with  cesium. 

This  report  describes  the  final  assembly  of  the  MPD  test  section,  the  design  and  construction 
of  the  fast  readout  unit  for  the  48- pin  electrodes,  and  the  results  of  the  runs  with  pure  helium. 

Experience  gained  from  the  runs  with  helium  seeded  with  cesium  is  discussed. 
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II.  RESUME  OF  PROGRESS 


MANUFACTURE  OF  CLOSED  LOOP  DEVICE 

During  this  report  period  the  diagnostic  closed  loop  MPD  device  was  completed. 

The  following  changes  were  made  in  the  system. 

•  The  charcoal  filters  were  removed.  Only  the  absolute  filters  remain  in  the  system. 

•  A  calibration  un-'t  for  the  heat  balance  flow  meter  was  designed,  installed,  and  checked 
out. 

FIRST  TEST  RUNS 

A  large  number  of  runs  with  pure  helium  were  made  during  which  pressure  probe  and  thermo¬ 
couple  measurements  were  obtained.  It  was  possible  to  establish  calibration  curves  for  all 
parameters  to  set  the  system  for  a  certain  mass  flow  at  a  given  temperature.  The  impurities 
were  cheeked  spectroscopically.  For  water  vapor  and  02.  continuous  sampling  techniques 
were  employed. 

To  date  four  runs  with  helium  seeded  with  cesium  have  been  made.  A  typical  but  preliminary 
value  for  a  gas  temperature  of  1000’K  was  1.0  volt  generated  voltage  and  0.4  ma  generated 
current. 

The  system  remained  operable.  More  runs  with  helium  seeded  with  cesium  are  scheduled.  It 
is  anticipated  that  after  completion  of  these  runs  a  few  changes  in  the  system  design  will  be 
made. 

THEORETICAL  INVESTIGATIONS 
Equilibrium 

A  new  equation  was  derived  describing  on  a  statistical  basis  the  ionization  equilibrium  of 
plasmas  in  magnetic  fields. 
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Nonoquilibrium 


Field  equations  and  transport  relations  for  a  reacting  multitemperature  plasma  are  presented 
in  this  report.  The  theory  involved  is  based  on  the  classical  approximation  neglecting  the 
perturbation  of  the  particle  distribution  function  by  the  reactions. 


4 


111.  THE  ALLISON  DIAGNOSTIC  CLOSED  LOOP  MPD  DEVICE 


FINAL  ASSEMBLY  OF  THE  MPD  TEST  SECTION 

Figure  1  is  the  schematic  diagram  of  the  MPD  test  section.  Figure  2  is  a  photographic  view 
of  the  section.  The  test  section  has  been  designed  specifically  for  diagnostic  investigation  of 
plasma  properties.  Emphasis  has  been  placed  on  versatility,  ease  of  maintenance,  and  sim¬ 
plicity. 


As  shown  in  Figure  1,  the  test  section  consists  of  a  flow  channel  electrode  block,  observation 
ports,  and  electromagnetic  assembly.  Figure  3  shows  the  tantalum  flow  channel  with  observa¬ 
tion  window.  The  plasma  enter  from  the  right.  The  vertical,  rectangular  pieces  are  the  elec¬ 
trode  ducts.  The  horizontal,  cylindrical  pieces  are  radiation  shields  which  cover  the  tantalum 
channel. 

The  tantalum  channel  and  the  electrode  ducts  are  lined  with  high  purity  alumina  panels  to  pro¬ 
vide  electrical  insulation  with  respect  to  currents  induced  in  the  7x  B  interaction.  They 
serve  also  as  protection  for  the  tantalum  walls.  The  linings,  shown  in  Figure  4,  are  easy  to 
replace.  Figure  4  is  a  view  into  one  of  the  electrode  ducts.  One  observation  window  can  be 
seen  f'-om  inside  the  duct.  The  rectangular  opening  which  is  also  visible  inside  is  the  exit  of 
the  plasma  channel. 

The  tantalum  duct  is  housed  in  a  stainless  steel  structure  containing  a  helium  protective  atmos¬ 
phere  to  prevent  oxidation.  Figure  5  is  a  close  view  of  the  MPD  Section.  Figure  6  is  an  assem 
biy  view  of  entire  section.  The  electrodes  consist  of  48  platinum  pins  shrunk  into  an  alumina 
block.  One  of  the  electrodes  is  shown  in  Figure  7.  The  two  pins  at  the  left  in  Figure  7  serve 
as  supports  for  a  Pt-PtRh  thermocouple.  An  electrical  lead  from  the  pin  is  connected  to  a 
corresponding  pin  of  a  50-pin  hermetically  sealed  connector.  The  connector  is  visible  in  Fig¬ 
ure  6. 

DIRECT  READOUT  UNIT 

A  direct  readout  unit  was  developed  to  display  the  generated  voltages  of  each  pin  separately  on 
an  oscilloscope  screen.  The  voltages  developed  between  each  pin  in  the  upper  electrode  and 
its  corresponding  pin  in  the  lower  electrode  are  indicated  by  the  oscilloscopic  display  depicted 

in  Figure  8. 
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Relatively  high  voltages  are  indicated  for  the  pin  pairs  of  the  No.  5  column.  Low  voltages  are 
indicated  for  the  right  half  of  the  first  pin-pair  row.  The  system,  as  presently  constructed, 
permits  scanning  of  all  voltages  in  approximately  one  second.  The  scanning  may  be  stopped 
and  the  voltages  appearing  across  the  pin-pairs  of  any  one  column  may  be  continually  displayed. 
The  circuitry  is  shown  in  Figure  9.  The  console  is  shown  in  Figure  10.  The  heart  of  the  scan¬ 
ning  unit  is  a  solenoid-driven,  24-position,  ten-pole  stepping  switch.  Alternate  positions  are 
grounded-- leaving  twelve  active  positions. 

Consider  row  No.  1  of  the  top  electrode.  The  twelve  pins  are  connected  by  cable  to  the  twelve 
active  positions  of  the  row  No.  1  Switch — Top.  The  twelve  pins  of  row  No.  1  of  the  bottom 
electrode  similarly  connect  to  the  active  positions  of  row  No.  1  Switch  —  Bottom.  Operation  of 
the  switch  causes  successive  pin  pairs  to  be  connected  to  the  differential  amplifier.  The  dif¬ 
ferential  amplifier  feeds  to  the  first  trace  of  the  oscilloscope.  A  Tektronix  oscilloscope  fitted 
with  a  four-input  plug-in  unit  is  employed.  The  remaining  rows  of  electrode  pins  provide  simi¬ 
lar  inputs  to  the  oscilloscope. 

It  is  desirable  to  place  variable  loads  across  the  pin  pairs.  These  are  provided  by  the  1000-ohm 
load  rheostats  as  indicated  in  Figure  9.  Knobs  for  adjusting  appear  in  Figure  10. 

It  is  also  desirable  to  be  able  to  connect  the  pins  in  various  series  and  parallel  combinations. 

To  facilitate  this,  the  ends  of  each  rheostat  are  brought  to  banana  plugs  on  the  front  panel. 

These  plugs  are  located  in  pairs  just  beneath  each  rheostat  adjusting  knob  as  indicated  in  Figure 

10. 

At  present,  the  stepping  switch  moves  at  only  one  speed — maximum.  The  jar  associated  with 
operation  of  the  solenoid  causes  contact  bounce  which  shows  up  as  noise  on  the  oscilloscope 
traces.  It  is  planned  to  provide  variable  rates  of  scan.  Slower  rates  of  scan  will  permit  the 
contact  bounce  to  die  out  before  the  next  step  is  taken. 

At  present,  scanning  is  initiated  by  depression  of  the  scanning  switch  shown  in  Figure  10.  This 
switch  applies  power  to  the  solenoid  and  also  delivers  a  trigger  signal  to  the  scope  for  initiat¬ 
ing  scanning  in  the  scope.  This  system  operates  very  well  for  a  single  rate  of  scan.  Introduc¬ 
tion  of  variable  rates  of  scan  makes  a  change  desirable.  It  is  planned  to  remove  the  scope 
triggering  signal  and  in  its  place  feed  the  scope  horizontal  sweep  from  voltages  derived  from 
the  stepping  switch.  This  will  ensure  that  each  voltage  appears  in  proper  position  regardless 
of  rate  of  scan. 
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Figure  1.  MPD  diagnostic  device. 
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Figure  1.  M PD  diagnostic  device. 


Figure  G.  Assembly  view  of  MPl)  section. 


Figure  7.  Electrode  consisting  of  48  platinum  pins 


12 


Note:  All  switches  ganged, 


Figure  9.  Circuitry  schematic. 
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Figure  10.  Direct  readout  unit. 
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IV.  TEST  RESULTS 


RUNS  WITH  PURE  HELIUM 

The  runs  with  pure  helium  served  as  a  shakedown  for  the  overaii  performance  of  the  system. 

The  main  point  of  interest  was  to  establish  calibration  curves  for  different  mass  flows  to  be 
abie  to  produce  a  given  temperature  by  varying  the  heater  power  input. 

Other  points  of  interest  were  the  performance  of  the  heater,  the  filters,  the  heiium  compres¬ 
sor,  and  the  cooling  system. 

A  iarge  number  of  runs  were  made  and  about  25  hours  of  running  time  were  accumulated. 
Between  the  runs  the  system  was  dismantled  for  inspection.  The  highest  temperature  was 
1500“K.  The  highest  mass  fiow  was  8.74  gm/sec. 

The  runs  with  pure  heiium  caused  no  attack  on  the  tantaium  parts  of  the  test  section.  Also, 
the  ceramic  parts  showed  no  attack.  However,  a  siight  deposit  on  the  ceramic  parts  was 
found.  Spectroscopic  analysis  made  of  the  deposit  indicated  W,  Ta.  Th.  and  traces  of  Sn  and 
Be.  Aii  these  materials  are  present  in  the  heater,  the  main  source  of  impurities.  Judging 
from  the  amount  of  deposit  and  the  running  time  involved,  one  can  estimate  that  these  impuri¬ 
ties  do  not  amount  to  more  than  1  ppm.  It  was  not  possible  to  detect  these  impurities  with 
spectroscopic  analysis  of  the  gas  during  the  runs. 

Measurements  of  water  vapor  content  and  oxygen  content  were  continuously  monitored.  The 
water  vapor  was  kept  down  to  10  ppm;  the  oxygen  content  varied  from  0.  5  to  i.  5  ppm. 

Static  and  total  pressure  probes  were  instaiied  to  obtain  calibration  curves  for  heater  power 
input  settings.  The  results  of  these  probe  measurements  are  shown  in  Figures  11  and  12. 

Figure  1 1  is  a  plot  of  gas  temperature  versus  heater  power  input,  while  Figure  12  gives  the 
gas  velocity  versus  heater  power  input.  Both  gas  properties  are  from  the  center  of  the  test 
section.  Figure  13  gives  a  piot  Mach  number  versus  power  input,  which  was  derived  from 
Figures  11  and  12.  it  is  intended  to  repeat  these  runs  with  thermocouples  in  piaee  of  the  probes 
and  compare  the  resuits. 
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•Temperature  measured  with  Thermocouple,  M  =  4.6 

Figure  11.  Gas  temperature  vs  heater  power  mass  flow  M  parameter 

1  iguro  N  shows  the  helium  jet  in  the  test  section.  The  electrodes  were  separated  sufficiently 
to  show  the  jet.  l'or  running  purposes,  electrode  distances  of  1  cm  and  2  cm  were  used. 

For  shakedown  purposes,  the  electrodes  were  connected  to  a  variable  load  resistor  and  output 
voltage  and  current  wore  measured.  For  runs  with  pure  helium,  of  course,  there  should  be 
no  output  at  all.  However,  a  small  output  voltatc  in  the  order  of  0.  2  v  and  a  small  output 
current  in  the  order  of  0.  3  ma  were  detected.  The  exact  values  depend  strongly  on  voltage  of 
the  heater,  mass  flow,  and  gas  temperature,  and  only  slightly  on  the  magnetic  field.  This 
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indicated  that  the  observed  output  was  due  to  a  leakage  current  from  the  heater  coming  either 
through  the  gas  or  along  walls.  If  it  is  a  leakage  current  through  the  gas.  its  value  may  be  in¬ 
creased  tremendously  when  injecting  cesium.  This  fact  has  to  be  borne  in  mind  when  one  is 
interpreting  data  from  power  runs  at  low  power  levels. 


Figure  12.  Gas  velocity  vs  heater  power 
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Figure  13.  Mach  number  vs  heater  power  mass  flow  M  parameter  (g/soc) 


A  more  extensive  investigation  of  tins  phenomenon  is  scheduled  and  the  results  will  be  in¬ 
cluded  in  the  next  Quarterly  Summary  report. 

HUNS  WITH  HELIUM  SEEDED  WITH  CESIUM 

The  objective  of  the  first  runs  with  the  seeded  working  fluid  was  to  learn  more  about  the  prob¬ 
lems  that  are  involved  in  the  operation  of  this  type  of  loop.  Major  points  of  concern  were: 

•  Material  attack  of  the  tantalum  structure 

•  Material  attack  of  the  ceramic  structure 

•  Possible  arcing  of  the  heater 

•  Performance  of  the  cesium  injection  system 

•  Performance  of  the  cesium  separation  and  seed  recovery 


Figure  14.  Helium  jet  in  test  section. 


Besides  this.  She  electrical  data,  even  ii  preliminary  and  crude,  give  tlx-  most  important  in¬ 
formation.  It  allows  the  design  of  tlie  proper  instrumentation  setup  and  defines  the  problems 
which  are  involved. 

The  preliminary  conclusions,  drawn  from  these  runs  were: 

•  The  system  is  still  operational;  therefore,  only  minor  design  modifications  arc  necessary 
due  to  the  experience  obtained  from  these  runs. 

•  Only  minor  arcing  did  occur  within  the  heater, 

•  \  judgment  on  material  attack  can  be  made  only  after  disassembly  of  the  system. 
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The  cesium  injection  system  has  to  be  improved.  Also,  the  use  of  a  two-stage  separator  in¬ 
stead  of  the  present  one-stage  separator  is  desirable.  The  quartz  window  devitrified  due  to 
the  cesium  attack  and  should  be  replaced  with  sapphire  windows.  These  changes  will  be  made 
as  soon  as  the  present  series  of  runs  is  completed. 

With  the  present  setup  the  generated  voltage  and  the  leakage  voltage  are  of  opposite  polarity. 
Therefore,  values  of  1.0  volt  and  0.4  ma  at  1000*K  gas  temperature,  reported  in  Seetion  11. 
are  pessimistic  values. 

\ftcr  learning  more  about  the  leakage  mechanism  and  after  more  runs  are  made  with  reversed 
magnetic  field,  it  will  be  possible  to  give  a  definite  value  of  the  generated  power  at  1000°K  gas 
temperature. 
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V.  TllEOKETiCAL  i  WEST  iG  AT  IONS 


iOMZATlON  EQUILIBRIUM  OF  PLASMA  ]\  MAGNETIC  FIELD 


Thermodynamic  Equilibrium  in  Plasma 


A  plasma  can  be  oniy  in  true  thermodynamic  equilibrium  if  it  is  homogeneous  and  if  escape  of 
radiation  is  negligible.  Consequently,  in  presence  of  a  magnetic  field  spatiai  constancy  of  this 
field  must  be  demanded,  inhomogeneity  or  existence  of  different  temperatures  in  the  com¬ 
ponents  is  necessarily  connected  with  irreversible  processes  like  diffusion,  heat  conduction, 
and  eiectricai  conduction,  which  disturb  the  thermodynamic  equilibrium.  It  is  obvious  that  the 
equilibrium  state  is  an  ideai  limit  to  which  an  actual  plasma  can  oniy  approximate. 

The  reactions  occurring  in  a  plasma  are  provided  by  inelastic  particie  coiiisons.  As  known 
from  statistical  mechanics,  the  contribution  of  these  interactions  to  the  totai  energy  can  be 
neglected  when  the  system  is  sufficiency  dilute,  in  the  Hamilton  operator  of  the  plasma  sys¬ 
tem, 


'Total 


.  ^  2 


j=i 


2  m. 


j 


where 


generalized  momentum  operator  of  the  j-th  particle 
charge  and  mass  of  the  j-th  particie 
interaction  energy  for  the  j-th  and  k-th  particle 
magnetic  moment  of  the  j-th  particle 


A  =  vector  potential  of  the  magnetic  field  B 


(1) 


therefore,  the  interaction  energies,  W,  of  the  particles  are  disregarded.  As  a  consequence, 
the  result  of  the  investigation  is  independent  of  the  iaws  of  interaction. 
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The  subject  of  the  investigation  is  a  homogeneous  plasma  submerged  into  a  constant  magnetic 
field  in  which  ionization  and  deionization  reactions  take  place  on  the  average  at  the  same  rate: 

a  *  atoms 

a  i  *  e,  e  =  electrons  (2) 

i  =  ions 


The  condition  for  thermodynamic  equilibrium,  formulated  in  terms  of  the  electrochemical 
potentials,  Mr,  of  the  components  r.  is 


r 


e 


M.  -M  *  0 

i  a  . 


(.1) 


With  7. r  as  the  sutne  of  states  and  \r  as  total  number  of  particles  of  the  r  -  component,  the 
potentials,  p  are  given  by 


-  k  T 


?ln  Zr 

ax,. 


T.  V.  Sr/ 


( I) 


The  sum  of  states  is  conveniently  expressed  by  the  collective  of  the  energy  eigenvalues,  K0. 
which  a  representative  particle  of  the  s-componcnt  can  exhibit:1 


■ 


V 


2  e'K«/k  T 


X,. 


This  relation  refers  to  weakly  interacting  particles  of  Bose  and  Fermi  types. 


Sum  of  States  for  Particles  in  Magnetic  Field 


(5) 


\ s  a  theoretical  model,  consider  a  plasma  box,  Figure  15,  of  side  lengths  Sx,  Sy,  S7  which 
contains  \0  electrons,  Xj  ions,  and  Xa  atoms  in  homogeneous  mixture.  The  magnetic  field, 
B,  is  assumed  to  be  directed  along  the  z-axis.  It  is  given  in  terms  of  the  vector  potential  by 


(6) 


Superscripts  denote  references  listed  in  Section  VI. 
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Figure  15.  Plasma  box  with  magnetic  field. 


The  plasma  is  treated  in  the  scmiclassical  approximation  which  is  based  on  the  assumption1 


Electron  Component 

The  magnetic  moment  of  the  electrons  due  to  the  spin  is  /T  =  cr  eh/m,,  where  e  >0 

2'  2 

is  the  elementary  charge.  As  the  spin  operator,  cr ,  commutes  with  the  Hamilton  operator  H 
and  the  constant  vector  potential,  \,  commutes  with  the  momentum  operator,  "p  the  wave 
equation  of  the  electrons  can  be  derived  ns: 
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The  wave  function  is  obviously  of  the  form 
i  (Px  x  t  Pz  z) 

,  h 

<t>  =  Const  X  c  Xu  (y) 


where  Px<  P  are  momentum  eigenvalues,  -  ®  <  Px,  Pv  <  ♦».  Separation  of  the  coordinates 
leads  to 


Equation  (9)  corresponds  to  the  wave  equation  for  a  linear  oscillator  of  energy,  Eogc  = 

.  e  n  pz"  n 

E  *  <r  n - - - .  oscillating  with  the  electron  gyration  frequency.  c  s  e  H/mc,  about  the 

me  2  m0 

point  yG  -P  /c  B.  Accordingly,  the  energy  eigenvalues  are’ 


•v  P  IT 
z' 


'  2/  m„  m„ 


2  mr 


’*  0,  1. 


(10) 


The  momentum  eigenvalues  Px  and  Pz  in  Equations  (9)  and  (10)  are  still  continuously  degener¬ 
ate.  This  degeneracy  is  readily  removed  upon  considering  the  finite  dimensions  of  the  box.  It 
results  for  the  number  of  momentum  eigenvalues  in  the  phase  space  element  APxl.x  -APZL,Z. 


At  = 


APx  K 

2  ”  h 


■IP  L 
z  z 

2  w  ti 


(11) 


Prom  the  condition  that  the  gyration  centers  must  be  contained  within  the  limits  0  iyi  Ly, 
namely  0  i  yQ  S  Ly,  it  follows  APX  *  e  BLy.  Consequently,  the  number  of  momentum  eigen¬ 
values  in  the  considered  phase  space  element  becomes* 


At 


V 

(2»rtl)2 


e  BA  Pz. 


V  5 


s  s  s 

x  y  z 


(12) 


By  extending  the  dimensions  of  the  box  to  infinity,  and  considering  Equations  (10)  and  (12),  it 
follows  for  the  electron  component: 
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Equivalent  to  this  condition  is  the  demand  for  the  collision  frequency  to  be  small  compared  to 
the  gyration  frequency  of  the  electrons, 


,coll« 

e  c 


(17) 


w  lie  re 


coll 


-  1 


k  T 

-  and«n  =  e  B/m„.  The  ease  of  strong  eollisional  perturbation 


of  the  electron  gyration  motion  is  not  regarded  here. 


Ion  and  Atom  Component 


Additional  energy  states  due  to  excitation  appear  in  the  ionic  and  atomic  component.  The  in¬ 
trinsic  energy  Es  of  a  particle  with  an  electronic  structure  in  state  s  is  determined  by  the  four 
quantum  numbers  ns,  Ls,  Ss,  Js,  when  external  fields  are  absent.  In  the  considered  magnetic 
field,  the  ion  or  atom  exhibits  in  each  state  s  due  to  the  quantized  orientation  of  its  magnetic 
moment,  /i(s)  =  Mgs  (eh/2  me) — in  a  statistical  sense— the  oreintational  energy  eigenvalues 
-  M  j»(s),  where  M  *  -  Js . t  Js.  The  nontranslational  energy  eigenvalues  of  the  con¬ 

sidered  particles  arc,  therefore. 


EM.  «  -  Es-*«s  * 


2  m. 


(18) 


It  is  assumed  here  that  the  spin-orbit  interaction  is  predominant.  The  Lande  splitting  factor 
is  defined  by 


.1  (,J  +  1 )  -»  (S„  ♦  1)  -  L  <L  +  1) 

S  S  S  o  5  o 

2,Is(Jg  ♦  1) 


( I  0) 


It  is  usual  to  write  Es  as  the  sum  of  the  energy  of  the  atom  in  ground  state,  «G.  and  the  ex¬ 
citation  energy  of  the  s-th  state,  k  s, 


E  =  «+  «,  K  sO  (20, 

.s  *o  “s'  “o 

By  consideration  of  Equations  (18),  (ID),  and  (20),  the  sum  of  state  and  chemical  potential  of 
the  ions  is  easily  derived  in  analogy  to  the  previous  case.  As  the  ion  gyration  radius  rj  = 

V  k  T  /c  B  is  Vm j/mc  times  larger  than  the  electron  gyration  radius,  two  situations  are 
of  interest: 
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1.  Unperturbed  Ion  Gyration  Motion 


Ion  gyration  radius  r;«  ion  mean  free  path!  j.  It  results 


A  -(> 

I  e  B  V  , - ^  \  2/  m,k 

1  —  - T~ 5  v2*m7kT  /  i  e  ‘ 
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'/  V*  V'  MRLt‘7JLVkT 
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M*  -Jg  8--0 


M‘  5  'kT  'n 


^Trrl^nEe  2'  m‘kT 

J  v  .n 


iL.,„  E  Ec'(-i"5"ih^'kT 


M*  -Jg  s=0 


2.  Perturbed  Ion  Gyration  Motion: 


Ion  gyration  radius  rj  >>  ion  mean  free  path/,.  It  results: 


,3/2  j 

( -ir~)  v  x  c  0 


M-  -  Jg  S  =  0 


M  ,  =  -  k  T  <  In 


1  m-  _  il 


M  =  -  s  =  0 
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Formally,  Equation  (23)  can  be  derived  by  carrying  through  the  limit  B  -*  0  in  the  translational 
part  of  Equation  (21). 

\n  expression  similar  to  Equation  (23)  is  valid  for  uncharged  partieles.  as  their  translational 
energy  states  are  unaffected  by  the  magnetic  field.  The  corresponding  expressions  for  the  sum 
of  states  and  chemical  potential  of  the  neutral  atom  component  are: 


The  double  sums  in  Equations  (21)  through  (2f>)  are  later  designated  by  ur  (B.  T)  (exritational 
partition  function  of  the  component  r,  r  -  i.  a). 

Ionization  Equations  for  Plasma  in  Magnetic  Field 


Tire  equilibrium  particle  composition  of  plasma  in  magnetic  field  is  readily  obtainable  by  in¬ 
serting  the  applying  chemical  potentials  into  the  condition  for  thermodynamic  equilibrium, 
Equation  (3).  Again,  the  following  eases  are  considered: 

1.  I'npcrturbed  Electron  and  Ion  Gyration: 

From  Equation  (3)  and  Equations  (15),  (22),  and  (20),  it  follows  as  ionization  equation: 


=  hJ 


^2  ir  rr>j 


( 2  n  ma) 


TFT 


_  /  l\  h  e  B 

[  e  B  ~|2  V  2/  mc  k  T 

L(2rrh)2J  ^ ° 

u  J  r  =  0 


E' 

v  =0 


JA  jl£  n 

2/mj  k  T 


y 


+  1/2 


-  a 

e 


tic  B 
me  k  T 


cr  =  -1/2 


»B(T) 

uB<T) 


(k  TV 


1/2 


(«o  -  «8> 


k  T 


(27) 


28 


29 


.Allison 


j)  hin 
M  -  0 


z 


lie  B 

-a - 


me  k  T 
e 


&  -  -  1/2 


2 


k)  lim 
B  -  0 


‘n 


(T) 


■  Is  (21>  * ')' 


I  k  T 


u«<T>. 


l.  a 


Thi*  limit  j  demonstrates  the  double  degeneracy  of  the  electrons  and  the  limit  k  the  (2Jg  ♦  11- 
times  degeneracy  of  the  ionic  and  atomic  electron  states  in  absence  of  the  field,  B  =  0. 

By  consideration  of  limits  l  through  k.  both  Equations  (27)  and  (28)  yield  in  the  limit  of  vanish¬ 
ing  field,  B  0: 


ne  nj 
na 


(30) 


This  is  the  well  known  Eggert-Saha  equation,  describing  the  ionization  equilibrium  of  field-free 
plasmas. 


It  is  immediately  relevant  that  the  ionizational  composition  of  plasma  in  magnetic  field  deviates 
from  the  Eggert-Saha  equilibrium  only  at  high  magnetic  field  intensities,  when  the  quanta  of  the 

c  B 

electron  gyration  motion.  tiw  =  h - ,  arc  of  the  order  of  magnitude  of  the  thermal  energy: 


k  T  ~ 


In  this  case,  the  magnetic  field  has  the  effect  of  remarkably  increasing  the  degree  of  ionization. 
The  condition  of  Equation  (31)  can  be  realized  by  means  of  high  intensity  short-time  electro¬ 
magnets  in  low  temperature  plasmas.  Under  these  extreme  conditions,  the  spin-orbit  inter¬ 
action  energy  in  the  atom  can  be  small  compared  to  the  energy  of  the  atomic  magnetic  moment 
in  the  exterior  field.  In  this  situation,  the  so-called  Paschen-Bach  effect  is  observed.  This 
phenomenon  influences  somewhat  the  form  of  the  excitational  partition  functions  as  defined  in 
Equations  (24)  and  (26). 
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It  can  be  presumed  that  the  quanta  of  the  ion  gyration  motion,  *  ti 


e  B 

- ,  are  small  compared 

mi 

to  the  thermal  energy,  kT,  at  all  technically  achievable  field  intensities.  Furthermore,  as 
the  influence  of  the  magnetic  field  is  of  interest  only  at  low  temperatures,  the  ions  and  atoms 
can  be  considered  to  be  in  their  ground  state.  It  follows  from  Equation  (27)  or  (28): 
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(32) 


where  the  hyperbolic  contangcnt  has  been  introduced.  Furthermore,  Ug  (Mc)  and  Up  (MQ)  dcsig 
nate  the  orientational  partition  function  of  the  atoms  and  ions  in  ground  state,  respectively. 
These  arc  obtained  from  the  corresponding  excitational  partition  functions  by  taking,  in  the  in¬ 
finite  scries  over  s,  only  the  first  term  (s  =  0).  When  Russell-Saundcrs  coupling  does  not 
apply,  the  partition  functions  have  to  be  modified  appropriately. 

Disregarding  nor.cquilibrium  effects,  the  ionization  degree  in  MPD  converter  plasmas  can  be 
calculated — in  a  very  rough  approximation — from  equilibrium  theory.  For  the  usual  region  of 
temperature  ( 1000-3000°K)  and  magnetic  field  intensity  (10-20  X  10^  gauss)  in  converter 
systems,  the  Eggert-Saha  equation  and  the  new  ionization  equation.  Equation  (32),  lead  to 
practically  the  same  result. 


Ql'ASl-KlNETIC  THEORY  OF  REACTING  PLASMA 

\  theoretical  description  of  multicomponent  converter  plasmas  should  Include  that  the  individual 
components  arc  not  in  thermal  equilibrium  with  each  other,  and  that  this  intcrcomponcnt  non¬ 
equilibrium  also  depends  strongly  on  the  inelastic  particle  collisons.  The  thermal  nonequilib- 
rhun  among  the  components  is  due  to  the  fact  that  because  of  the  large  mass  differences  of  the 
different  particles,  approximate  equilibrium  in  the  single  component  is  much  more  easily 
achieved  than  between  the  light  and  heavy  particle  components.  The  possible  nonequilibrium 
to  be  observed  in  stationary  state  is  determined  primarily  by  the  mechanism  of  electrical 
energy  input,  the  distributive  intercomponcnt  convective  processes,  and  the  reactive  particle 
flux  and  reactive  energy  flux  sources. 
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\  kinetic  theory  for  multitemperature  plasmas  based  on  the  Boltzmann  equation  has  been  given 
by  Sehradc*  and  was  later  reconsidered  in  more  detail  by  Schindler*.  This  theory  is  extended 
in  the  following  paragraphs  to  reacting  converter  plasmas  in  a  classical  approach  negleciing 
the  perturbation  of  the  particle  distribution  functions  by  the  reactions.  In  this  connection  it  is 
remarkable  that  in  an  actual  plasma  these  perturbations  are  presumably  smoothed  in  the  region 
above  the  thermal  velocity  by  microinstabilities. 


Theoretical  Method 

Designate  by  ws  *  v*g  -  v*s  the  actual  particle  velocity  of  the  s-component  as  observed  in  a  sys¬ 
tem  moving  with  the  mean  mass  velocity.  "vs.  of  the  s-component.  In  the  individual  frame  of 
reference  of  the  s-component.  the  n-tli  order  velocity  moment  of  the  s-component  is  then  de¬ 
fined  by 

t)  •  dVg  (33) 

The  first  thirteen  scalar  velocity  moments,  which  have  a  simple  physical  meaning,  arc: 
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Scalar  pressure,  ps,  and  thermal  energy  density,  Us,  arc  defined  by  the  reduced  second  order 
velocity  moment: 


r  n  _  2 

ps  =  msJw;fsdvs  =  -  Us 

By  multiplication  of  Boltzmann's  Equation  with  the  tensor  mc  wl”*- ;  =  m„  w„  ,•  w=  . .  and 

“  3i  1J  >  •  •  S  S  ,  1  S(  J 

subsequent  integration  with  respect  to  velocity  space,  it  can  be  shown  that  the  n-th  order 
velocity  moments  satisfy  the  general  equation  of  change  (summation  according  to  Einstein  con¬ 
vention), 
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The  expressions  on  the  right  side  of  Equation  (34)  represent  the  collision  integrals  for  elastic 
and  reactive  collisons,  respectively. 

The  elastic  particle  collisions  are  treated  exclusively  as  two-body  collisions.  For  elastic  two- 
bodv  interactions,  the  collisional  change  of  the  distribution  function  is  given  by  7 
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The  reactive  many-body  collision  integrals  arc  treated  in  a  phenomenological  way  by  introduc¬ 
ing  reaction  velocities. 

Neglecting  the  influence  of  the  reactive  collisions,  the  velocity  distribution  functions  of  the 
components  s  are  given  in  the  so-called  1  3- moment-approximation  by* 
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From  Equation  ( 34).  the  field  and  transport  equations  for  the  s-eomponents  are  obtained  as 
momentum  conservation  equations. 

The  elastic  collision  integrals  in  the  momentum  conservation  equations  can  be  readily  evaluated 
by  means  of  Equation  (3G)  for  the  following  models: 

•  Interaction  corresponding  to  rigid  elastic  spheres  of  radii  r8  and  rr  for  encountei  s  be¬ 
tween  a  neutral  and  another  neutral  particle,  or  between  a  charged  and  a  neutral  particle. 

In  the  center  of  mass  system  of  the  colliding  particles,  the  differential  cross  seetion  is 


a 


ar 


(37) 


In  evaluating  the  collision  integrals.  Equation  (37)  leads  to  a  total  transport  cross  section 

i 

(Hnmsauer  cross  section) 


QSr  '  '(rs  +  rr)2 


(38) 


•  Interactions  corresponding  to  Coulomb  scattering  fields  for  encounters  between  particles 
of  charge  es  and  er,  where  the  collision  parameter  is  cut  off  at  the  Debye  shielding  dis- 
tance, 


In  the  center  of  mass  system  of  the  colliding  particles,  the  differential  cross  section  is 
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In  evaluating  the  collision  integrals,  Equation  (39)  leads  to  a  total  mean  transport  cross  section 
(Grosdovcr  cross  section)7 

rr  T  °s  er  "1  ^  ^  ^  sr  ®"1  , 

Qsr  ■  1  [r^rj  '"'-■'s'-1  1+  l^~\ 

The  reduced  mass,  /isr,  and  reduced  temperature,  Tsr.  are  defined  by 


s'  ms  +  mr 


T  =  M 
1 sr  rsr 
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\s  reactions,  ionization  and  recombination  processes  are  considered.  It  is  assumed  that  the 
ionization  of  the  neutral  atoms  is  provided  by  electron  collisions,  and  the  recombination  of  the 
ions  is  due  to  three-body  electron- ion-electron  collisions.  The  corresponding  reaction  cqun- 


c  +  a  ; — *  i  +  e  +  c 

This  model  applies,  in  particular,  to  the  situation  prevailing  in  low  temperature  (T  5  3000°K) 
alkali  seeded  (seeding  ratio*  ~  10'2)  noble  gas  plasmas  of  atmospheric  pressure.  In  speci¬ 
fying  this  type  of  plasma,  the  following  indices  arc  introduced: 

Index  o;  for  noble  gas  atoms 

Index  a:  for  neutral  seeding  atoms 

Index  i:  for  ionic  seeding  atoms 

Index  c:  for  electrons 

In  the  considered  temperature  range,  the  noble  gas  component  is  treated  as  nonreacting. 

The  zero  order  reactive  collision  integrals  represent  the  reactive  mass  production  of  the 
species  s.  For  the  regarded  reactions,  introduce  r+  as  ionization  reaction  velocity  and  r  as 
recombination  reaction  velocity.  According  to  the  definition  of  reaction  velocities,  the  follow¬ 
ing  relations  can  be  established: 


/•.  -reactive 
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it  is  for  the  electron  component  s  =  e:  (r,  P)  -  (a,  o),  (i.  e);  and  appropriately  for  tire  re¬ 
maining  components.  The  *s  designate  the  stoichiometric  coefficients  in  the  reaction  equation 
as  formulated  in  Equation  (42): 


=  *2-1  =  tl.  =  »l.  k3  =  -t.  »-0  =  0  (44) 

It  is  immediately  recognized  that  summation  of  Equation  (43)  with  respect  to  the  s-componettts 
yields  zero  (conservation  of  mass  in  the  reactions). 


The  total  reaction  velocity,  f  =  r*  -  r  *.  is  given  in  terms  of  the  (electron  eoiiision)  ionization 
probabiiitv  coefficient.  S0.j,  and  (three-body  eieetron-ion-eieetron)  recombination  probability 
coefficient,  Uejc.  by 
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(45) 


The  reaction  coefficients  thus  defined  are— in  first  approximation— functions  of  the  electron 

» 

temperature  aione.  According  to  Eiwcrt  : 


(46) 


“e.e  =  r*’ 
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where  the  correction  function  of  magnitude  of  order  I  is  given  by 

The  symbols  in  Equations  (45)  and  (46)  have  the  meaning: 


(48) 
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e  =  2.998  X1 01 0 


=  6.  6  X  iO’25  crn^  =  Thomson  scattering  cross  section 

=  Sommcrfeld  fine  structure  constant 
cm 

-  =  velocity  of  light  in  free  space 

sec 
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(  n  =  number  of  electrons  in  the  highest  shell  n  of  the  atom  in  ground  state 

*  n  *  ionization  energy  of  the  atoms  from  n-th  shell 

*  II  1  13.  53  e  v  -  ionization  energy  of  hydrogen  atom  in  ground  state 

h2  _8 

:l  s  ,  2  2  =  529  x  1°  cm  *  first  Bohr  radius  of  hydrogen  atom 

1*  mp  e“ 

5  cxcitational  portion  function  of  the  ions  (uj=  2  jJ,  ♦  1  at  low  temperatures) 
un  =  cxcitational  partition  function  of  the  atoms  (ua  =  2  j“  ♦  1  at  low  temperatures) 
f  =  factor  taking  into  account  ionization  from  excited  states 

It  is  remarkable  that  Equations  (45)  and  (4G)  are  derived  under  the  assumption  of  a  Maxwellian 
electron  distribution,  and  that  the  recombination  reactions  are  considered  to  be  in  equilibrium 
with  the  ionization  reactions.  The  use  of  equilibrium  reaction  coefficients  in  the  description 
of  substantial  nonequilibrium  plasmas  has  to  be  understood  in  the  sense  of  an  approximation. 

As  to  the  reactive  collision  integrals  for  the  first  order  velocity  moments,  it  is  assumed  that 
their  contribution  to  the  momentum  transport  in  the  components  is  negligible.  The  reactive 
collision  integrals  for  the  reduced  velocity  moments  of  second  order,  representing  the  reactive 
energy  production  in  the  s-component,  arc  obviously  of  the  form: 


It  designates  AES  the  reactive  variation  of  energy  of  the  s-particles.  In  the  considered  re¬ 
al  lion  mechanisms,  the  ionization  energy  is  taken  primarily  from  the  electron  component,  the 
recombination  energy  is  fed  primarily  into  the  electron  component.  Furthermore,  with  every 

3 

transformed  s-particle,  the  energy  vg  -  k  Tg  is  liberated  in  the  mean.  Consequently, 


According  to  the  present  state  of  theory,  reaction  velocities  for  nonequilibrium  plasmas  arc 
not  available. 
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AE0  =  0 

Summation  of  Equation  (49)  with  respect  to  the  s-components  yields  only  zero  for  T  =  r  .  t 
in  the  trivial  case  of  thermodynamic  equilibrium.  In  nonequilibrium,  therefore,  the  plasma 
as  a  whole  exhibits  a  reactive  energy  source. 


Field  Equations 


The  field  equations  for  the  components  are  obtained  by  evaluating  the  general  equation  of 
change.  Equation  (34),  for  the  zero  order  velocity  moment,  the  first  order  velocity  moment 
and  the  reduced  second  order  velocity  moment.  The  appearing  reactive  collision  integrals  are 
given  in  Equations  (43)  and  (49).  After  the  explicit  calculation  of  the  elastic  collision  integrals  , 
the  results  can  be  summarized  as  follows  (s,  r  *  e,  i.  a,  o): 


Equations  of  Mass  Conservation: 
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Equations  of  Momentum  Conservation: 

+  V  •  *ps  -  ns  es  (E  ♦  7.  Xlll  * 


d  z.  I  =. 

-  *  v-  ‘  T  vs 

at  J 


sr  (vr  -  vs>  * 


r/s 

/J  V  o  Msr3/2  (  ?r-  - 

2sa sr  ns  nr  sryYT^r\nrmr  ns  ms/ 


(51) 


(52) 
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Equations  of  Energy  Conservation: 
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The  single  coefficients  depend  on  the  law  of  interaction  and  arc,  explicitly,  for: 
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Higid  Sphere  Interactions: 


Coulomb  Interactions: 


The  physical  meaning  of  the  expressions  on  the  right  side  of  Equations  (51),  (52),  and  (53)  is  in 
most  cases,  from  their  structures  self- relevant.  In  Equation  (51),  the  terms  proportional  to 
the  total  reaction  velocity  T  =  T+  -  r  represent  the  reactive  mass  production  of  species  s  per 
cm1  and  per  sec.  The  intcrcomponent  convective  terms  in  Equation  (52)  arc  the  friction  and 
thermodiffusion  forces  between  the  different  components.  In  Equation  (53),  the  term  propor¬ 
tional  to  the  total  reaction  velocity  T  =  T  +  -  r  '  is  the  reactive  energy  source  in  the  s-oom- 
ponent.  The  terms  proportional  to  the  temperature  difference  ATps  =  Tr  "  Ts  represent  the 
heat  transfer  between  the  components  of  unequal  temperature.  The  following  term  represents 
the  thermalization  of  the  kinetic  energy  of  the  interpenetrating  components  due  to  the  friction 
forces.  The  last  term  considers  the  power  of  the  thermodiffusion  forces. 
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Transport  Relations 


The  field  equations  are  to  be  supplemented  by  the  transport  equations  for  the  anisotropic  part 
of  the  pressure  tensor,  *s  =  Pg  -  pg  U.  and  heat  current  vector  q^. 


Pressure  Tensor 


1  lie  general  transport  equation  for  the  anisotropic  part  of  the  pressure  tensor  is  obtained  bv 
applying  the  general  equation  of  change.  Equation  (34).  to  the  corresponding  velocity  moment. 
1  o  make  this  transport  equation  convenient  for  applications,  terms  associated  with  heat  cur¬ 
rents  and  reaction  influences  and  terms  nonlinear  in  velocity  gradients  are  neglected.  Also, 
specializing  for  phenomena,  which  are  slow  compared  to  the  relaxation  time,  rg,  for  viscous 
momentum  transfer,  the  final  relation  is: 


-  2  «  v  X  i 
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(54) 


Note  that  the  cyclotron  frequency  vector,  <-g  =  es  B/ms,  depends  on  the  sign  of  the  respective 
charge,  c8,  The  relaxation  times  are  for.  it'1 


ltlgid  Sphere  Interactions: 


Qss  from  Equation  (38); 
Coulomb  Interactions: 


Qss  from  Equation  (40). 


The  vector  product  between  a  tensor  and  a  vector  is  defined  by 

(TXA),j  *  Tki‘jklAi 
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where  <jkl  is  the  permutation  tensor.  For  two  vectors  A  and  B,  the  symmetric  traceless 
tensor  is  defined  as 


A  B>>iJ  =  !  (Ai  V  Bi  Aj  '  3  Sij  Ak  Bk) 


As  an  example,  consider  a  plasma  with  the  magnetic  field  in  the  y-direction.  In  this  case,  the 
solution  of  Equation  (55)  is: 
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=  Vyx  ’  ‘  ,  **2r2  [Us.^-Us,yZ“sTs] 


2  2 

zz  +  us,  xx^  “*  “s  Ts  +  us,  zx  2  "s  Ts 


rr  -it 

s.xy  s 


JT  =  ir 

s,  xz  s.zx 


ir  =17 

s.yz  s.zy 


=  -  - >~T~;  k  xz  +  \  (us,  xx  '  US,  zz>  2  -S  Ts| 

1  +  -»«sT8  L  J 

:  '  9  2  [Us-  >z  +  Us’  xy  “s  Ts  I 

l  +“iTs  L  j 


The  viscosity  has  been  defined  in  the  usual  way. 


**8  5  1’S  TS 


The  tensor  components  us  ^  are  abbreviations  for  the  components  of  the  tensor  «t7vs». 
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Heat  Current 


The  general  transport  equation  for  the  heat  current  vector  is  obtained  by  applying  the  general 
equation  of  change.  Equation  (34),  to  the  corresponding  velocity  moment.  As  in  the  previous 
case,  various  simplifications  are  introduced.  Neglected  are:  (I)  terms  proportional  to  the 
gradients  of  the  velocity  fields:  (2)  terms  connected  with  gradients  of  the  nonisotropic  pressure 
tensor:  and  (3)  reaction  influences.  Restricting  again  to  phenomena,  which  are  slow  compared 


to  the  characteristic  relaxation  time  rg  = 


3 

2 


rg,  Equations  (55-50 


it  follows  for  the  regarded 


transport  equation: 


q  -  q  X  *  r  * 
S  S 


5  k  *  ,  T 

-  -  -  Ps  %  VTS  *  r  S( 

2  ms  8  8  s  s 


s  s 


(GI) 


The  analytical  solution  of  this  equation  is  explicitly: 


«s‘-I.w2  r‘2{  (As  VTs  *  Ts  ®s)  *  [.  rs  '  (Xs  VTs  '  Ts  ,)],  Ts 

-s%^(xs  vT.-':n) 


(62) 


where  the  thermal  conductivity  is  defined  as 


5  k 

*8  E  2  rn^  Tf 


(63) 


The  forces  which  produce  the  thermal  currents  arc  temperature  gradients,  VTS,  and  diffusion 
forces.  {  .  of  interpenetrating  components  with  different  particle  masses  (Dufour  effect)* 


h  --  12  ns  nr 


UrkTsr) 


3/2 


*sr  ^vr  '  vs> 


(64) 


r^s 


According  to  the  applying  law  of  interaction,  the  coefficients  are:* 


Rigid  Sphere  Interactions: 

8  _Tj_  \  jl_ 

*sr  1  T  9  T  ^ 

J  1 sr  c  1  sr 
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Coulomb  Interactions: 


The  thermodiffusion  effect,  liquation  (6-4),  is  of  the  intercomponent  convective  type  and.  there¬ 
fore.  of  primary  importance  for  the  possible  nonequilibrium  which  can  establish  between  the 
components. 

Further  Remarks 


In  the  presented  description  of  reacting  multitempcraturc  plasma,  excitation  reactions  are  in¬ 
cluded  only  by  means  of  a  corrective  factor.  In  principle,  a  more  exact  treatment  of  excitation 
reactions  and.  consequently,  also  of  ionization  reactions  from  excited  levels,  is  possible  as 
the  reaction  coefficients.  Equations  (4G)  and  (47),  can  be  applied  to  any  level.  However,  the 
number  of  equations  would  be  increased  considerably,  as  alone  for  every  particle  kind  in  ex¬ 
cited  state  n  (n  =  1,  2,  3,  ..  .to  )  a  separate  conservation  equation  has  to  be  formulated.  For 
this  reason,  following  Elwcrt*  a  phenomenological  factor.  £  .  correcting  for  the  excitation  re¬ 
actions  has  been  introduced.  In  low  temperature  plasmas,  £  —*  1  can  be  assumed  when  the 
atoms  are  practically  in  the  ground  state. 


Another  point  is  the  convergence  of  the  scries  development  of  the  velocity  distribution  function. 
Equation  (3G).  Obviously,  the  anisotropic  pressure  components  must  be  small  compared  to  the 
static  pressure, 


«  Ps  • 


and  the  heat  current  vector  components  must  satisfy  the  condition. 


ms 


1 


3/2 


In  continuouslv  operated  converter  plasmas,  these  relations  arc,  in  general,  satisfied,  as  the 
gradients  of  the  velocity  and  temperature  fields  can  be  considered  to  be  small.  In  a  shock  wave 
structure,  with  strong  inhomogoncitics  of  the  fields,  for  comparison,  the  conditions  may  well 
be  violated.  The  perturbation  of  the  distribution  functions  by  the  reactions,  finally,  should  be 
negligible  in  the  considered  temperature  range. 

Reference  to  this  part  of  the  report  will  be  made  in  various  applications  which  arc  in  progress. 
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